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Observation of charge asymmetry dependence of pion elliptic flow and the possible 

chiral magnetic wave in heavy-ion collisions 
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We present measurements of 7r“ and tt^ elliptic flow, V 2 , at midrapidity in Au+Au collisions at 
= 200, 62.4, 39, 27, 19.6, 11.5 and 7.7 GeV, as a function of event-by-event charge asymmetry, 
Ach, based on data from the STAR experiment at RHIC. We find that 7r“ (tt^) elliptic flow linearly 
increases (decreases) with charge asymmetry for most centrality bins at = 27 GeV and higher. 

“ 200 GeV, the slope of the difference of V 2 between 7r“ and tt"*" as a function of Ach 
exhibits a centrality dependence, which is qualitatively similar to calculations that incorporate a 
chiral magnetic wave effect. Similar centrality dependence is also observed at lower energies. 

PACS numbers: 25.75.Ld 


In heavy-ion collisions at the Relativistic Heavy Ion 
Collider (RHIC) and the Large Hadron Collider (LHC), 
energetic spectator protons produce a strong magnetic 
field reaching eBy « [T], or ^ 3 x 10^"^ T. The in¬ 


terplay between the magnetic field and the quark-gluon 
matter created in these collisions might result in two phe¬ 
nomena: the chiral magnetic effect (CME) and the chiral 
separation effect (CSE). The CME is the phenomenon 
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of electric charge separation along the axis of the mag¬ 
netic field in the presence of a finite axial chemical po¬ 
tential [IHS]. The STAR Hi] and PHENIX (TOjlll] Col¬ 
laborations at the RHIC and the ALICE Collaboration 
at the LHC [H] have reported experimental observations 
of charge separation fluctuations, possibly providing ev¬ 
idence for the CME. This interpretation is still under 
discussion (see e.g. and references therein). The 

CSE refers to the separation of chiral charge, which char¬ 
acterizes left/right handedness, along the axis of the mag¬ 
netic field in the presence of the finite density of electric 
charge [mm- In this Letter, we report the results from 
a search for these effects using a new approach. 

In a chirally symmetric phase, the CME and CSE 
can form a collective excitation, the chiral magnetic wave 
(CMW). It is a propagation of chiral charge density in 
a long wavelength hydrodynamic mode [T5I120) . The 
CMW, which requires chiral symmetry restoration, man¬ 
ifests itself in a finite electric quadrupole moment of 
the collision system, where the “poles” (“equator”) of 
the collision system acquire additional positive (negative) 
charge [TH|. This effect, if present, will increase(decrease) 
the elliptic flow of negative (positive) particles. Elliptic 
flow refers to an azimuthally anisotropic collective motion 
of soft (low momentum) particles. It is characterized by 
a second-order harmonic in a particle’s azimuthal distri¬ 
bution, (/, with respect to the reaction plane azimuthal 
angle, 'Jrp, which is determined by the impact parameter 
and the beam direction, 

V2 = (C0S[2((() - ifRp)]). (I) 

The CMW is theoretically expected to modify the elliptic 
flow of charged particles, e.g. pions, on top of the baseline 

^base(7r±) p] 

P2(^±)=u^^‘’(7r±)T^Aeh, (2) 

where r is the quadrupole moment normalized by the 

net charge density and = (A^+ — N_)/{N^ + N_) 
is the charge asymmetry of the collision system. As the 
colliding nuclei are positively charged, the average charge 
asymmetry, (Ach), is always positive. Thus, the Ach- 
integrated V 2 of 7r“ (7r+) should be above (below) the 
baseline because of the CMW. However, the may 

be different between 7r+ and Tr~ because of several other 
possible physical mechanisms [HHHj- It is preferable to 
study CMW via the Ach dependence of the pion V 2 other 
than Ach-integrated V 2 - 

This Letter reports the Ach-differential measurements 
of the pion V 2 , based on Au-|-Au samples of 2 x 10® events 
at 200 GeV from RHIC year 2010, 6 x 10^ at 62.4 GeV 
(2010), 10® at 39 GeV (2010), 4.6x 10^ at 27 GeV (2011), 
2 X 10^ at 19.6 GeV (2011), 1 x 10^ for 11.5 (2010) and 
4 X 10® for 7.7 GeV (2010). All events were obtained 
with a minimum-bias trigger which selects all particle- 
producing collisions, regardless of the extent of overlap 



FIG. 1. (Color online) (a) the distribution of observed charge 
asymmetry from STAR data and, (b) the relationship between 
the observed charge asymmetry and the charge asymmetry 
from HIJING generated events, for 30-40% central Au-I-Au 
collisions at 200 GeV. In this centrality, the mean charge 
asymmetry (Ach) of HIJING events is about 0.004. The errors 
are statistical only. 

of the incident nuclei [5^. Charged particle tracks with 
pseudorapidity I??] < 1 were reconstructed in the STAR 
Time Projection Chamber (TPC) [5^. The number of 
charged particles within I 77 I < 0.5 is used to define the 
centrality. The centrality definitions and track quality 
cuts are the same as those used in Ref. HZj, unless oth¬ 
erwise specified. Only events within 40 cm (50 cm for 
11.5 GeV and 70 cm for 7.7 GeV) of the center of the de¬ 
tector center along the beam line direction are selected. 
To suppress events from collisions with the beam pipe 
(radius = 3.95 cm), a cut on the radial position of the 
reconstructed primary vertex within 2 cm was applied. A 
cut on the distance of the closest approach to the primary 
vertex (DCA < 1 cm) was applied to all tracks to sup¬ 
press contributions from weak decays and/or secondary 
interactions. 

The observed Ach was determined from the measured 
charged particles with transverse momentum px > 0.15 
GeV/c and \r]\ < 1; protons and anti-protons with 
Pt < 0.4 GeV/c were excluded to reject background pro¬ 
tons from the nuclear interactions of pions with inner 
detector materials. Fig. [^a) shows an example of the 
observed Ach distribution, which was divided into five 
samples roughly containing equal numbers of events, as 
indicated by the dashed lines. Fig. [^b) shows the rela¬ 
tionship between the observed Ach and the Ach from the 
HIJING event generator [5B|, where the same cuts as used 
in data were applied to calculate Ach- The relationship 
is linear. To select pions with high purity, we eliminate 
charged particles more than 2cr away from the expected 
energy loss of pions in the TPC. For energies less than 
or equal to 62.4 GeV, elliptic flow measurements were 
carried out with the V 2 {t] sub} approach |29j . where two 
subevent planes register charged particles with rj > 0.3 
and r] < —0.3, respectively. Pions at positive (negative) 
r] are then correlated with the subevent plane at nega¬ 
tive (positive) 77 to calculate V 2 - The 77 gap of 0.3 unit 
suppresses several short-range correlations such as the 
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FIG. 2. (Color online) (a) pion V2{2} as a function of ob¬ 
served charge asymmetry and, (b) V 2 difference between 7r“ 
and TT^ as a function of charge asymmetry with the tracking 
efficiency correction, for 30-40% central Au-|-Au collisions at 
200 GeV. The errors are statistical only. 


Bose-Einstein interference and the Coulomb final-state 
interactions m- There are correlations that are unre¬ 
lated to the reaction plane that are not suppressed by 
the r] gap, e.g. those due to back-to-back jets. These 
are largely canceled in the V 2 difference between 7r“ and 
7r+. For 200 GeV, the two-particle cumulant method 
V2{2} [301 EH was employed, which was consistent with 
V 2 {ri sub}, and allowed the comparison with the n2{4} 
method discussed later in this letter. The same rj gap 
was also used in the V2{2} analysis. To focus on the soft 
physics regime, only pions with 0.15 < pt < 0.5 GeV/c 
were used to calculate the pT-integrated V2 , and this pr 
range covers 65-70% of all the produced pions. The cal¬ 
culation of the pT-integrated V2 was corrected with the 
p^-dependent tracking efficiency for pions. 

Taking Au-|-Au 200 GeV collisions in the 30-40% cen¬ 
trality range as an example, the pion V 2 is shown as a 
function of the observed Ach in Fig. |^a). The 7r“ V 2 
increases with increasing observed Ach while the tt"'" V 2 
decreases with a similar magnitude of the slope. After 
applying the tracking efficiency to Ach, the V 2 difference 
between 7r“ and has been fitted with a straight line as 
shown in Fig. j^b). The slope parameter, r, from Eq. 2, 
is positive and qualitatively consistent with the expecta¬ 
tions of the GMW picture. The fit function is non-zero 
at the average charge asymmetry (Ach), which is a small 
positive number in case of Au -I- Au collisions. This indi¬ 
cates the Ach—integrated V 2 for 7r“ and 7r+ are different, 
which was observed in Ref. [33] . We follow the same pro¬ 
cedure as above to extract the slope parameter, r, for 
all centrality bins at 200 GeV. The results are shown 
in Fig. 1^ together with simulations using the UrQMD 
event generator [53] and with the theoretical calculations 
with GMW [Mj with different duration times of the mag¬ 
netic field. For most data points, the slopes are positive 
and reach a maximum in mid-central/mid-peripheral col¬ 
lisions, a feature also seen in the theoretical calculations 
of the GMW. The gray bands in Fig. include three 
types of systematic errors: the DGA cut for pion tracks 
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FIG. 3. (Golor online) The slope parameter, r, as a function of 
centrality for An-I-An collisions at 200 GeV. Also shown is the 
UrQMD |33| simulation, and the calculations with GMW |34| 
with different duration times. The grey bands include the 
systematic errors due to the DGA cut, the tracking efficiency 
and the pr range of particles involved in the event plane deter¬ 
mination. The cross-hatched band indicates the STAR mea¬ 
surement with the U2{4} method and the height of this band 
shows only the statistical error. 


was tightened to 0.5 cm, to study the contribution from 
weak decays, which dominates the systematic errors; the 
tracking efficiency for charged particles was varied by rel¬ 
ative 5%, to determine the uncertainty of Ach; and the 
Pt range of particles involved in the event plane determi¬ 
nation was shrunk from [0.15, 2] GeV/c to [0.7, 2] GeV/c, 
to further suppress short-range correlations. The Ach bin 
center may not accurately reflect the true center of each 
Ach bin in Fig. as the V 2 measurements are effectively 
weighted by the number of particles of interest. Such an 
uncertainty on r has been estimated to be negligible for 
most centrality bins, except for the most peripheral col¬ 
lisions, where this systematic error is still much smaller 
than the statistical error. 

To further study the charge-dependent contribution 
from jets and/or resonance decays, we separated positive 
and negative particles in each subevent to form positively 
(negatively) charged subevents. Then each 7r+ (tt”) is 
only correlated with the positive (negative) subevent in 
the opposite hemisphere. The slope parameters thus ob¬ 
tained are statistically consistent with the previous re¬ 
sults though with larger uncertainties. 

The event plane reconstructed with particles recorded 
in the TPG approximates the participant plane, the mea¬ 
sured V 2 are not the mean values, but closer to the root- 
mean-square values |35j . Another method, f2{4} |36j is 
supposed to better represent the measurement with re¬ 
spect to the reaction plane. For 20 — 50% Au-|-Au col¬ 
lisions at 200 GeV, the slope parameter obtained with 
V2{4:} is illustrated with the cross-hatched band in Fig.[^ 
which is systematically lower than the f 2 { 2 } results, but 
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FIG. 4. (Color online) The slope parameter r as a function 
of centrality for all the collision energies under study. For 
comparison, we also show the calculations with CMW [34] 
with different duration times. The grey bands carry the same 
meaning as those in Fig. 


still has a finite positive value with a larger statistical 
error. 

Since the prediction of the consequence of CMW on 
V 2 [IHl mil this subject has recently drawn increasing 
attention from theorists [SU [37ti42] . It was pointed out 
in Ref. (42] that local charge conservation at freeze-out, 
when convoluted with the characteristic shape of V 2 {pt) 
and ^ 2 ( 77 ), may provide a qualitative explanation for the 
finite V 2 slope we observe. Such an effect depends on the 
strength of the Ach dependence on the mean pT and the 
77 -dependence of V 2 - However, our measurements were 
carried out in a narrow px range ([0.15,0.5] GeV/c) and 
with a {pT){Ach) variation of 0.1% at most. Furthermore, 
the measured 77 -dependence of V 2 is only half as strong 
as that used in Ref. [12]. We estimate the contribution 
of this mechanism to be smaller than the measurement 
by an order of magnitude. 

To check if the observed slope parameters come from 
conventional physics, such as Coulomb interactions, or 
from a bias due to the analysis approach, we carried out 
the same analysis in Monte Carlo events from UrQMD. 
As shown in Fig. [^ the slopes extracted from UrQMD 
events of 200 GeV Au-I-Au collisions are consistent with 


zero for 10-60% centrality collisions, where the signal is 
prominent in the data. Similarly, the AMPT event gener¬ 
ator mi |44| also produces events with slopes r consistent 
with zero. With the AMPT model, we also studied the 
weak decay contribution to the slope, which was negli¬ 
gible. On the other hand, the CMW calculations [T5] 
demonstrate a similar centrality dependence of the slope 
parameter. Recently, a more realistic implementation of 
the CMW mi suggested that the CMW contribution to 
r is sizable, and the centrality dependence of r is sim¬ 
ilar to the data. In these theoretical calculations such 
centrality dependence mainly results from the central¬ 
ity dependence of the magnetic field and the system vol¬ 
ume. Quantitative comparisons between data and theory 
require further work on both sides to match the kine¬ 
matic regions used in the analyses. For example, the 
measured Ach only represents the charge asymmetry of 
a slice (I 77 I < 1 ) of an event, instead of that of the whole 
collision system. We expect these two values of Ach to 
be proportional to each other, but the determination of 
the ratio will be model dependent. In addition to the 
UrQMD and AMPT simulation studies which reveal no 
trivial correlation between Ach and pion V 2 , tests were 
performed using the experimental data. For example, 
Ach and the pion V 2 were calculated in two kinematically 
separated regions, i.e., different rapidity bins. In such 
cases, the slope parameters decrease but remain signifi¬ 
cant and positive. This may reflect the local nature of 
the Ach dependence of V 2 , but additional theoretical de¬ 
velopment is necessary. 

Figure [^ shows a similar trend in the centrality de¬ 
pendence of the slope parameter for all the beam energies 
except 11.5 and 7.7 GeV, where the slopes are consistent 
with zero with large statistical uncertainties. It was ar¬ 
gued [2T| that at lower beam energies the Ach-integrated 
V 2 difference between particles and anti-particles can be 
explained by the effect of quark transport from the pro¬ 
jectile nucleons to mid-rapidity, assuming that the V 2 
of transported quarks is larger than that of produced 
ones. The same model, however, when used to study 
772 (tt”) — U 2 ( 7 r+) as a function of Ach, suggested a nega¬ 
tive slope US], which is contradicted by the data. 

The mean field potentials from the hadronic 
phase [22] and the partonic phase [24] also qualitatively 
explain the Ach-integrated V 2 difference between parti¬ 
cles and anti-particles, especially at lower beam energies. 
In general, the mean field potential is expected to be 
positively correlated with Ach and thus may explain the 
trends in those data, but no conclusive statement can be 
made here due to the lack of specific predictions. This 
effect may be tested in the future by studying the V 2 
slopes, whose V 2 ordering is opposite to that of 7 r=*=. 

In summary, pion V 2 exhibits a linear dependence on 
Ach, with positive (negative) slopes for 7 r“ ( 7 r+). The 
772 (tt”) — 772 ( 7 r+) increases as a function of Ach, qual¬ 
itatively reproducing the expectation from the CMW 
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model. The slope r of W 2 (^ch) difference between 7r“ 
and 7r+ has been studied as a function of centrality, 
and we observe a dependence also similar to the calcu¬ 
lation based on the CMW model. The slope parameter 
r remains significantly positive for 10 — 60% centrality 
Au-I-Au collisions at = 27 — 200 GeV, and displays 
no obvious trend of the beam energy dependence with 
the current statistics. None of the conventional models 
discussed, as currently implemented, can explain our ob¬ 
servations. 
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